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Time-reversal symmetry is important to optics. Optical processes can run in a forward or backward 
direction through time when such symmetry is preserved. In linear optics, a time-reversed process of laser 
emission can enable total absorption of coherent light fields inside an optical cavity of loss by time-reversing 
the original gain medium. Nonlinearity, however, can often destroy such symmetry in nonlinear optics, 
making it difficult to study time-reversal symmetry with nonlinear optical wave mixings. Here we 
demonstrate time-reversed wave mixings for optical second harmonic generation (SHG) and optical 
parametric amplification (OPA) by exploring this well-known but underappreciated symmetry in nonlinear 
optics. This allows us to observe the annihilation of coherent beams. Our study offers new avenues for 
flexible control in nonlinear optics and has potential applications in efficient wavelength conversion, 
all-optical computing. 

Time-reversal symmetry dominates in many physical systems; it allows a physical process to reverse in 
a backward direction of time. This powerful symmetry enables many practical applications under such 
reversible principle, with examples including acoustic and electromagnetic focusing using time- reversal 
mirror 1 " 4 , phase conjugation mirror in optics 5,6 , etc. Especially in optics, light transmission in random medium 
can be enhanced by several orders by exploring the time- reversal symmetry 7 " 9 , similarly for enhancing second 
harmonic signals in random medium 10 . Recently, a novel concept named "coherent perfect absorber" (CPA) 11 " 13 
that explores time-reversed process of laser emission has shown that incident coherent optical fields can be 
perfectly absorbed by a time-reversed optical cavity, where the gain is replaced with an equal amount of loss. Also, 
the incident fields and frequency should coincide with those of corresponding lasing modes with gain under time 
symmetry. Recent efforts have studied CPA properties with different geometries 11 " 18 , though most in the linear 
optics regime. A nonlinear version of CPA has been theoretically proposed to investigate signal and idler beams' 
phase-varying dynamics in the presence of a pumping beam under a time-reversed optical parametric oscillation 
(OPO) scheme 14 . These studies of time symmetry have been attracting increasing attention, since they provide 
alternative and substantial ways to manipulate light in the nonlinear regime. 

In this work, we experimentally exam the time-reversal symmetry for two classical nonlinear wave-mixing 
processes-SHG and OPA, characterize their nonlinear properties as opposite to their time-reversal counterparts, 
and reveal the nontrivial dynamics of phase varying in time-reversed nonlinear wave-mixing schemes. These 
backward-nonlinear wave mixings lead to the unique property of annihilation of coherent beams in a nonlinear 
quadratic medium by time reversal. Unlike the case of CPA in the linear regime where incident fields are totally 
absorbed and converted into heat 1114 , annihilation of incident fields can lead to the generation of new fields. Such 
backward-parametric interactions may have a future application in efficient wavelength conversion for better 
long-wavelength detection, e.g. mid-IR, THz. More interestingly, a flexible phase control can be utilized to probe 
nonlinear dynamics during wave mixing, and redirect wave mixing forward or backward in time. This offers new 
techniques for flexible control in nonlinear optics and has potential applications in all-optical computing. 

Results 

SHG and time-reversed SHG. Figure 1 shows the scheme for time-reversed SHG. Two identical thin beta barium 
borate crystals (denoted BBOl and BB02 hereafter) cut for Type I SHG@1064 nm are used. The crystal length L is 
short enough (100 um) to ensure pump non-depletion (small signal approximation), i.e. A x (z) = A^O). The corres- 
ponding phase-mismatching vectors Ak (A/c') in BBOl (BB02) can be independently tuned by rotating the crystals. The 
intensity of the second harmonic (SH) after BBOl is 21 : 1 2 = 2d 2 ^a)lll/(nin 2 SoC 2 )L 2 smc 2 (AkL/2) = I norm smc 2 (AkL/2\ 
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Figure 1 | The schematic of SH generation and annihilation. The 

wave-mixing behavior of SHG process in the first nonlinear medium is 
time-reversed in the second (outlined in dashed box). The condition for 
time-reversal symmetry is obtained by adding a n phase shift to FW/SH 
phase difference and by nipping the sign of phase-mismatching vector ( Ak' 
= — Ak). For a spontaneously grown SHG (starting from zero) in BBOl, its 
time-reversed process leads to "perfect" SH annihilation in BB02 under 
time-reversal symmetry. 

where I is the wave intensity. We refer the normalized SH intensity 
Inorm to unity hereafter, which equals the output of BBOl at the 
phase-matching condition, rif is refractive index and e 0 is me 
vacuum permittivity. In BBOl, SH is generated as AkL/2 is 
controlled in the range of (— n, n), and the phase difference 
between the generated SH and FW is Acp = AkL/2. Thus BBOl is 
both an SH source generator and a phase controller. Both the pump 
FW and generated SH are then incident into the second nonlinear 
medium- BB02 after passing through a thin quartz plate for phase 
shift. For spontaneously-grown SHG (starting from zero) in BBOl, 
the time- reversed process points in exactly the opposite direction: SH 
signals with the appropriate relative amplitude and phase to the 
pumping beam can be totally annihilated when they incident onto 
a nonlinear crystal. 

We set BBOl for several different conditions where AkL/2 = 
— 0.37571, —0.22571, 0, 0.3571, 0.67L Figure 2 shows the total SH output 
after BB02 measured by scanning the phase mismatching vector Ak\ 
Significant annihilation of SH in BB02 occurs when time-reversal 
symmetry requirements are met, i.e. when the dips of the SH output 
curves are located at Ak' = —Ak. These results clearly indicate that, 
for a low conversion efficiency SHG, there always exists a symmet- 
rically reversed SHG process that is capable of "perfect absorbing" its 
SH wave when coherently illuminated by both FW and SH. The 
corresponding attenuation of each dip of the curves is measured to 
be 11.5 dB, 12.4 dB, 14.5 dB, 11.4 dB and 10.6 dB, and yet the SH 
generated in BBOl can, in theory, be totally absorbed in BB02 11,14 . 
Such modulation depth is limited by experimental conditions, e.g. 
reflection, pumping intensity, wavefront, coherence 13 . Moreover, 
complete annihilation cannot be achieved if quantum noise is 
considered 1424 . 

OPA and time-reversed OPA. OPA can be regarded as a reversed 
wave mixing to SHG from a quantum picture: OPA splits one photon 
into two while SHG combines two photons. However, they are not 
inherently time reversal to each other, for example, the initial 
pumping beams differ. Compared with reversed SHG, the 
pumping effect can play a critical role in reversed OPA geometry. 
Here we extend our study to a reversed two-channel OPA, which 
qualitatively requires two coherent input beams and a pumping 
beam. Annihilation of waves can only be achieved when the 
relative phase, amplitude and pumping level conditions are 
reached, as schematically shown in Fig. 3a. 
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Figure 2 | Experimental demonstration of reversed SHG. The measured 
total SH intensity with different FW/SH inputs into BB02 produced by 
varying its phase-mismatching vector Ak\ The dips are the signature of 
CPA, where SH generated in BBOl is significantly absorbed by time- 
reversed SHG in BB02. (Insets: Photographs of the SH spot at indicated 
conditions.). 

In Fig. 3b and 3c, we show the theoretical and experimentally 
measured parametric gain when the phase difference between the 
three waves is varied. Far from Ag? = n/2, the system behaves as an 
OPA device (<I> > 1). Note that at A(p = n/2, a dip with parametric 
gain O < 1 is observed in the gain curves, indicating Signal/Idler 
annihilation, when pumping above the threshold of OPA. Such a dip 
is the clear signature of "colored" OPA-CPA 14 , which is also the time- 
reversed process of its counterpart OPA above threshold. In an ideal 
case, the product of the gain factor for the in-phase and out-phase 
signal would be a unity. This is a straightforward outcome of time- 
symmetry rule. However, in our experiment, there are plenty of 
conditions that lead us to this result. One critical parameter is the 
pulse duration. The actual pulse of a pump wave with shorter time 
duration does not completely overlap with Signal/Idler pulses, which 
is the major reason for discrepancy between the experimental and 
theoretical curves, especially around the lowest attenuation points. 
Nevertheless, these dips still clearly indicate that a reversed-OPA 
device is capable of attenuating Signal/Idler given coherent illumina- 
tion with correct phases 26 , which behaves as a "colored" OPA-CPA 
system 14 . The result demonstrates two operating states of a wave- 
mixing process in its forward and time-reversed direction. 
Nevertheless, in our experiment, the gain depletion regime is never 
reached. In this regime, OPA-CPA is still possible when energy can 
be drained completely from the pumping beam 14 ' 27 . Here further 
study is required. 

Discussion 

The SHG process converts a fundamental wave (FW) into its second 
harmonic (SH) through the nonlinear response of a medium. How- 
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Figure 3 | The schematic of OP A and its reversal, (a) The schematic of Signal/Idler pair amplification and annihilation, (b) Behavior of the parametric 
gain O vs. total phase difference A(p for various pump intensities, (c) The experimentally measured overall O vs. A(p. Gain variety along pump 
pulses is taken into account in data fitting. 



ever, the energy flow can oscillate between the fundamental and 
harmonic waves, determined by the phase mismatch condition and 
phase difference of the interacting waves. In an undepleted-pump 
scheme, we consider the quasimonochromatic waves with carrier 
frequencies of fundamental wave at co l5 and second harmonic wave 
at <x> 2 = 2cd x . After representing the electric field as E { (z,i) = 
Ei(z) exp ( — iojjt) + ex. = A,-(z) exp[ — i(cOjt — kjz)] + ex., where i = 
1, 2 refers to FW and SH, respectively, and where A,-(z) is the slowly 
varying amplitude (slow- varying envelope approximation), the wave 
equation can be expressed as 21 : 



dA 2 
dz 



ioj\d e j 
k 2 c 2 



k 7 and d e jj is the 



where the phase-mismatching vector is Ak = 2k x — * 2 
effective nonlinearity, and c is the speed of light. Generally, in the 
SHG process SH signals grow along the propagation in a nonlinear 
optical crystal. However, they may also decay due to the familiar 
problem of phase-matching. Here the phase plays an important role: 
it directs the time axis where energy flows. For example, for SHG 
schemes in quasi -phase-matching (QPM) gratings, one can purposely 
flip the 2nd- order susceptibility (x (2) ) of ferroelectric domains thus 
adding a phase jump to SH waves in order to prevent SH falling back 
to FW due to dephasing. This encourages SHG conversion. 
Fundamentally, this originates from an interference effect similar to 
the linear CPA, but such interference is nonlinear one due to non- 
linear polarization between SH and FW beams at molecular level 19 . 
Thus, phase difference is crucial when considering the time-reversed 
process of SHG as well. 

To consider the time-reversed counterpart, we take the complex 
conjugate of Eq. 1 to get 



dA* 2 
dz 



k 2 c 2 



-A? 



JAkfz 



(2) 



Note that Eq. 2 is the same as Eq. 1 except for an extra phase differ- 
ence (A(p = cp 2 — 2<pi = n) and a sign-flipped phase-mismatching 



vector (Ak' = —Ak). The complex conjugate fields E*(z), which 
represent the backward fields, satisfy Eq. 2. With such exact time- 
reversal configurations, SHG beams undergo a backward process 
with respect to the counterpart in the forward-time direction under 
a time reversible environment that excludes magnetism and loss. For 
a non- seeded scheme, SH starts to grow from quantum fluctuations 
draining the energy from FW along the propagation. The exact time 
reversal of the problem should lead SH to zero, as it is the initial stage 
for SHG, similar to CPA in the linear case. Here the model is based on 
a monochromatic wave or waves with simple envelopes. For more 
complex waves, one has to time-reverse waveforms in space and time 
as well, e.g. through time-reversed mirrors 22 or nonlinear wave 
mixings 23 . 

To gain insight into the nonlinear wave-mixing properties of the 
proposed time- reversed structure, we calculate the total output after 
BB02 by integrating Eq. 1 along the propagating path. The total SH 
output is studied by scanning phase-matching vectors for both crys- 
tals with/ without a n phase shift to Acp, as shown in Fig. 4. The points 
of interest are where the SH generated in BBOl is totally cancelled 
out in BB02 leaving only FW out of BB02. In Fig. 4a, where two 
crystals are simply cascaded (without phase shift), SH only vanishes 
if AkL/2 = Ak'L/2 = ±n/2 (The blank points at the four corners are 
not the desired ones, and this is due to non-SH input and the absence 
of SH generated by BB02.). It is well-known that SHG at phase- 
mismatching condition induces energy flow oscillation between 
FW and SH at the period of twice of the coherent length L c = nl 
Ak, one can also purposely manipulate SHG output by phase con- 
trol 28 . However, in order to achieve the total absorption of SHG, we 
have to consider time-reversed scheme. Figure 4b shows the typical 
SH intensity along the propagating path at the conditions of phase- 
matched, quasi-phase-matched, and phase-mismatched conditions 
(red, green and blue lines, respectively). Note that SH only flows back 
under the phase mismatching condition, which occurs twice over the 
coherent length. And along its propagation, SH exhibits symmetry 
between the two crystals. This stimulates the quest to explore 
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Figure 4 | The total SH output characteristics and the SH intensity along the propagation path, (a) The normalized total SH output without a % phase 
shift between the two crystals. The two white dots show two situations when SH generated in BBOl is cancelled out in BB02 due to dephasing. (b) SH 
intensity along the propagation length in the two crystals in phase-matching and mismatching situations, (c) The normalized total SH output is 
greatly modified after introducing an additional n phase shift to A(p in between the two crystals. The diagonal dashed line shows total SH outputs equal to 
none, indicating any SH generated in the first nonlinear medium can always be cancelled out in the second, (d) SH intensity along the propagation length 
in the two crystals. The SH intensity evolution manifests itself in spatially symmetric pattern, which is also an indication of time-reversal of the wave 
mixing process. White dash dots and line in (a) and (c) indicate zero intensity. 



time-symmetry to re-converting SH back to FW at any moment 
despite their current phase relations. 

According to the analysis above, time-reversal symmetry of SHG 
can be implemented by introducing an additional n phase difference 
before BB02 and reversing the phase mismatching vector in BB02. 
We re-plot the SH output after BB02 in Fig. 4c. The SH generated in 
BBOl can always be cancelled out in BB02 whenever A A;' = — Ak is 
fulfilled. This reveals a condition where SHG generated by BBOl 
undergoes annihilation along a time-reversed path in BB02. 
Figure 4d illustrates this idea better by showing SH intensity along 
the propagating path across both crystals. With these conditions, 
the intensity profiles of SH preserves the spatial symmetry along 
the propagation axis through both crystals. Even more convincing 
is that this symmetry is preserved under the phase-matched 
condition in Fig. 4d. It is also worth to mention that the n phase 
jump is introduced purposely to trigger time-reversal symmetry 



using a tilt thin quartz plate sandwiched between two BBO crystals. 
In a nonlinear manner, such phase shift can arise nonlinearly in a 
"cascading nonlinearity" 29 , where an effective kerr-like x (2) (3oo;2oo,co) 
nonlinearity is created through cascading two x (2) (2co;oo,co) pro- 
cesses inside one single second-order nonlinear crystal. However, 
in this geometry, time symmetry is not preserved; second harmonic 
beams are converted into third harmonics instead of the fun- 
damental ones like in our case here. In all these circumstances, the 
phase plays a core role in determining the direction of the wave- 
mixing process. 

Then we consider non- depleted Type II degenerate OP A of three 
quasimonochromatic waves with carrier frequencies of fundamental 
wave at co p , signal and idler at w s = co { = co p /2. By slow- varying 
envelope approximation, neglecting group velocity mismatch 
(GVM) and group velocity dispersion (GVD), the envelops satisfy 
the coupled wave equations 21 
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(3) 



(4) 



The subscripts p, s, i stand for pump, signal and idler, respectively. 
The phase mismatching vector Ak = k p — k s — k t is required to be 
zero for an efficient OPA device. One can find the phase-matched 
parametric gain of signal/idler with respect to pump intensity and 
total phase difference (Aq? = 2(p p — q? s — (p t ) to be 14 ' 21,25 : <!> = 
\cosh (gz) — /exp( — iA(p) sinh(gz)| 2 , in which g = 2^/w s Wi/(n s rii) 
d e ffAp/c is proportional to the pump amplitude. The results with 
different pump intensities are plotted in Fig. 2b. Gain and attenu- 
ation are governed by the total relative phase Aq>. Hence, the relative 
phase matters most when we consider the time-reversed OPA. Since 
Ak = 0, OPA and its time-reversed process is accessible in the same 
condition, and an OPA device can behave like an OPA-CPA device 14 . 
Theoretically, the in-phase signal can be amplified by the factor O, 
whereas the out-of-phase signal is attenuated by the same factor. The 
dip in the gain curve at A(p = nil is a clear indicator of annihilation of 
Signal/Idler beams when pumping above the threshold of OPA. 

In conclusion, we have demonstrated time -reversed wave mixings 
for SHG and OPA. This enables us to observe the annihilation of 
coherent beams under time-reversal symmetry by varying the rela- 
tive phase of the incident fields. Time-reversed SHG is able to absorb 
coherent wave at second harmonic frequency of its pump. For time- 
reversed OPA, we show that the OPA can simultaneously amplify 
and attenuate coherent signal and idler waves when pumped above 
threshold. Our study provides a versatile platform for flexible control 
in nonlinear optics and potential applications in efficient wavelength 
conversion, all -optical computing. 

Methods 

Experiment setup for reversed SHG. The laser source used is an all-fiber 
femtosecond oscillator delivering 250 fs pulses centered at 1064 nm at a repetition 
rate of 80 MHz. A long-focal-length lens is used to loosely focus the light beam. And 
the two BBO crystals are placed on both side of the focus where the wavefront is also 
symmetrical. The desired n phase shift is obtained by inserting a tilt thin quartz plate 
in between. The tilt angle is calculated according to the dispersion relation such that 
SH wave accumulates an additional n phase shift than FW when passing through the 
plate. 

Experiment setup for reversed OPA. The gain media consist of a pair of anti-parallel 
BBO crystals pumped by an e-polarized frequency- doubled beam (532 nm) from a 
Q-switched nanosecond laser at a repetition rate of 20 Hz. The spatial wall-off of the 
three waves accumulated in the first crystal will be compensated by the second anti- 
parallel one. The crystals are cut for wavelength-degenerate type II OPA scheme 
pumped at 532 nm and are anti-reflection coated at both 532 and 1064 nm. The 
fundamental wave delivered by the laser is nominal 4 ns in duration (FWHM). And 
the pump beam, which is its second harmonic at 532 nm, has a nominal FWHM 
duration of 3 ns. A telescope system is used to reduce the pump beam diameter to 
2 mm. The fundamental wave (1064 nm) is attenuated to the extant that is much 
weaker than the pump (< 1/100), and its beam diameter is also reduced by a telescope 
system to a slightly smaller diameter to overlap with the pump. A half-wave plate is 
used to rotate the FW beam by 45 degrees to make the orthogonally polarized signal 
and idler input equal in amplitude on the onset of the crystal. 
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